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Abstract 

The  elastic-scattering  of  electrons  from  atomic  helium 
in  the  ground  state  is  investigated.  It  is  shown  that  for  low 
energy  incident  electrons  the  scattering  problem  reduces  to  solving 
an  ordinary  integro-di fferential  equation  for  the  scattering 
wave -function.  A  method  is  discussed  to  obtain  approximate  solutions 
to  the  integro-differential  equation  by  variational  principles. 

The  extremum  condition  of  the  variational  method  is  formulated  into 
a  general  N  x  N  matrix  equation  which  reduces  to  a  2  x  2  eigen-value 
matrix  problem  for  the  phase-shift  of  the  scattering  electron. 

An  algorithm  is  presented  to  obtain  the  collislonal  cross-section 
for  elastic  scattering  as  a  function  of  incident  electron  energy^ 


I 


Introduction 


Background 

It  Is  well-known  that  approximate  methods  based  on  variational 
principles  are  among  the  most  powerful  techniques  used  in  solving 
quantum  mechanical  problems.  In  bound-state  calculations,  the 
Raleigh-P.itz  variational  method  gives  energy  eigen-values  which  are 
correct  to  second  order  in  the  error  of  an  assumed  trial  wave-function. 

The  method  also  determines  the  free  parameters  of  the  wave-function 
and  guarantees  that  the  calculated  eigen-values  are  upper  bounds 
to  the  exact  energy  eigen-values  of  the  bound  system. 

In  scattering  theory  however,  the  total  energy  of  the  system 
is  specified  in  advance.  Variational  principles  are  used  to  determine 
to  first  order  the  free  parameters  of  the  trial  wave-function  which 
contain  all  information  of  the  scattering  event.  Unfortunately,  no 
general  bounding  conditions  are  yet  known  on  the  solutions  of  the 
continuum  problem  calculated  by  variational  methods.  Thus  there  exists 
no  defined  standard  (e.g.,  upper  or  lower  bounds)  with  which  to  compare 
computed  results  using  different  trial  wave- functions. 

In  spite  of  this  difficulty,  variational  principles  applied  to 
scattering  probUr-  offer  several  advantages  over  other  approximation 
methods.  First,  variational  methods  require,  in  general,  far  less 
numerical  effort  than  non -variational  methods.  Second,  results 
obtained  by  variational  calculations,  although  not  internally  bounded, 
compare  favorably  with  those  obtained  by  more  powerful  non-variational 
techniques.  Finally,  physically  complex  interactions  such  as  polarization 
and  electron  correlation  effects  can  easily  be  included  in  a  variational 
calculation  within  the  assumed  forms  of  the  trial  wave -funct ions . 

It  is  for  these  reasons  that  variational  principles  are  still 
competitive  with  other  more  powerful  non-variational  scattering 
techniques . 

The  first  serious  attempt  to  apply  variational  principles 
to  electron-helium  scattering  was  made  by  Moiseiwitsch  (Ref.  1). 
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Using  Che  variational  method  proposed  by  Hulthen  (Ref.  2)  and  a 
two  parameter  trial  wave -function,  Moiseiwitsch  found  good  agreement 
(107.)  for  the  S-wave  elastic  scattering  phase  shift  with  the  known 
numerical  solutions  of  the  same  integro-differential  equations  performed 
earlier  by  Morse  and  Allis  (Ref.  3).  Moiseiwitsch  further  suggested 
that  improved  results  could  be  obtained  by  more  complex  trial 
wave - f unc  t ions . 

In  a  detailed  study  of  variational  principles  applied  to 
continuum  problems,  Schwartz  (Ref.  4)  first  noted  that  at  certain 
system  energies  variational  calculations  possess  intrinsic  anomalous 
behavior.  He  showed  that  as  more  complex  trial  wave -functions  and 
more  adjustable  parameters  are  added  to  the  variational  calculations, 
the  solutions  do  not  converge  uniformly,  and  can  turn  out  totally 
erroneous  for  specific  incident  electron  energies.  Schwartz  also 
studied  the  nature  of  these  anomalies.  He  showed  that  although  a 
trial  set  of  functions  is  linear  with  respect  to  a  set  of  corresponding 
parameters,  the  quantity  which  is  extremi*ed  is  the  expectation  value 
of  the  scattering  Hamiltonian  operator.  This  quantity  is  always  quadratic 
in  the  parameters  of  the  trial  wave -function.  As  a  result,  the 
parameters  of  the  wave-function  have,  in  general,  two  roots,  only  one 
of  which  is  physical.  At  certain  system  energies  these  roots  can 
be  complex  and  the  variational  method  will  break  down. 

Several  attempts  have  been  made  to  eliminate  the  anomalous 
nature  of  a  single  variational  method  by  combinations  of  two  different 
variational  methods.  Conditional  calculations  proceed  in  the  same 
manner  as  a  single  variational  method.  However,  when  anomalous  results 
begin  to  appear  in  one  variational  calculation,  a  transfer  is  made  to 
a  second  variational  method.  Although  two  independent  methods 
possess  intrinsic  energy  anomalies,  they  do  not,  in  general,  occur 
in  the  same  energy  regions.  Thus  by  alternating  back  and  forth 
whenever  a  single  method  fails,  accurate  results  have  been  obtained 
over  the  entire  spectrum  applicable  to  the  physical  limitations  of 


the  approximations  made. 

Malik  (Ref.  5)  and  Rudge  (Ref.  6)  proposed  combining  the 
variational  method  of  Hulthen  and  Kohn  (Ref.  7).  Application  of  this 
technique  has  been  only  applied  to  electron-hydrogen  scattering  by 
Rudge  (Ref.  6).  Nesbet  (Ref.  8)  proposed  the  method  of  combining  the 
variational  principles  of  Kohn  and  the  Inverse  Kohn,  Rubinow  (Ref.  9). 
Both  combinational  methods  demonstrated  the  anomaly-free  nature  of  the 
technique  but  were  restricted  in  practice  to  elastic  scattering  for 
hydrogen,  since  the  methods  could  not  be  generalized  to  include 
inelastic  scattering. 

The  most  successful  anomaly-free  variational  technique  is 
the  Transformation  Method,  Harris  and  Michels  (Ref.  10).  This 
method  introduces  a  transformation  which  ensures  that  the  anomalous 
energy  regions  are  forced  outside  the  range  of  the  energy  spectrum 
of  interest.  Application  of  the  transformation  technique  to  variational 
calculations  of  electron  elastic  scattering  from  helium  was  first 
performed  by  Michels,  Harris  and  Scolsky  (Ref.  11).  This  work  was 
extended  by  Sinfailam  and  Nesbet  (Ref.  12).  Trial  exponential 
wave -functions  with  decreasing  numerical  indices  were  used.  The  results 
of  the  calculations  of  the  cross-sections  compared  extremely  well 
(3%  estimated  probable  error)  to  the  experimental  measurements  of 
Golden  and  Bandel  (Ref.  13). 

Nesbet  (Ref.  14)  refined  the  calculations  of  Sinfailam  and 
Nesbet  (Ref.  12)  for  electron  elastic  scattering  from  helium.  The 
results  for  the  cross-section  compared  extremely  well  (1.57.)  with  the 
best  non-variational  technique,  the  R-matrix  Method,  O'Malley,  Burke 
and  Berringtor.  (Ref.  15).  The  Nesbet  (Ref.  14)  calculation  represents 
to  date  the  best  variational  results  on  the  electron-helium  elastic 
scattering  croos -section. 

The  single  drawback  of  the  anomily-free  variational  methods 
is  the  increase  in  numerical  labor  required  above  the  simpler  variational 
calculations.  Recently,  Shankland  (Ref.  16)  proposed  a  variational 
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method  which  retains  the  simplicity  of  the  single  variational  approach. 
Recognizing  that  a  single  variational  method  fails  only  in  small 
energy  regions  near  non-physical  anomalous  poles,  he  proposed  a 
simple  interpolating  technique  to  obtain  an  accurate  cross-section 
in  the  singular  regions.  If  successful,  this  method  will  reduce 
the  numerical  labor  required  for  computing  accurate  elastic 
electron-atomic  cross-sections  over  a  large  energy  spectrum. 

Objective 

In  the  present  work,  the  Shankland  Variational  Method  is 
formulated  for  the  electron-helium  elastic  scattering  problem. 

Specifically,  an  algorithm  is  developed  to  obtain  variational 
solutions  for  the  scattering  phase-shifts  as  a  function  of  incident 
electron  energy. 

Approach 

In  Chapter  II  the  basic  formulation  of  the  electron-helium 
scattering  problem  is  presented.  In  particular,  the  close-coupled 
formulation  is  discussed.  Specifically,  the  simplest  case  of  the 
closed-coupling  method,  the  one-state  approximation,  is  developed 
in  detail.  In  Chapter  III  the  Matrix-Variational  Method  is  formulated. 
First  the  general  variational  method  is  applied  to  the  electron-hel ium 
scattering  equation.  Second  the  results  of  the  variational  formulation 
are  developed  into  a  matrix  eigen-value  problem.  In  Chapter  IV  a 
detailed  algorithm  is  presented  for  the  calculation  of  the  elastic 
scattering  phase-shifts  as  a  function  of  incident  electron  energy. 

Finally,  an  interpolation  method  is  discussed  for  the  special  case 

in  which  the  Matrix-Variational  Method  mathematically  yields  anomalous 

results.  Appendix  A  is  a  compilation  of  the  A-Matrix  elements 

required  in  the  computation.  Appendix  B  is  a  compilation  of  the 

B-Matrix  elements  also  required  in  the  computation.  Appendix  C 

is  a  compilation  of  the  definite  integrals  used  throughout  the  calculation. 
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II. 


The  Quantum  Mechanical  Problem 


Basic  Formulation 

In  this  section  the  quantum  mechanical  formulation  of  electron 
impact  with  helium  will  be  presented.  The  basic  framework  for  this 
approach  can  be  found  in  a  number  of  articles,  Massey  (Ref.  17), 
Moiseiwitsch  and  Smith  (Ref.  18),  as  well  as  texts  by  Mott  and 
Massey  (Ref.  19)  and  Geltmann  (Ref.  20).  In  the  following,  the 
formulation  presented  by  Mott  and  Massey  (Ref.  19)  is  the  most 
convenient.  Atomic  units  will  be  used  throughout,  (Ref.  19). 

The  wave  equation  for  the  scattering  of  an  electron 
by  a  helium  atom  is  given  as 

(2.t; 

\  r,  r,  r,  rlx  ^  > 


where  are  the  distances  of  the  three  electrons  from  the 

helium  nucleus ,  rl3 ,  are  the  distances  between  the  electrons, 
and  £  is  the  total  energy  of  the  system.  The  total  system 
wave -function  is  an  anti -symmetric  function  in 

the  spin-space  coordinates  of  each  pair  of  electrons. 

The  standard  method  of  attacking  :  he  general  N'-bodv 
scattering  problem  is  to  expand  the  total  svstem  wave-function 
in  terms  of  the  orthonormal  set  of  target  wave- funct ions .  For 
electron-helium  scattering  we  simply  have 

oo 

YCtofaxAvy)-  I  A  ^(r,<r,,fi<rak)Fn(r3<ra)  (2.2) 

nso 

where  denotes  the  n^  state  system  wave-function  for  the 

unperturbed  helium  atom.  It  is  an  anti-symmetric  wave-function 
in  the  spin-space  coordinates  of  the  helium  electrons  and  is  a 
solution  to  the  helium  equation 
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f  v«a  4  ?x*  a  +  it  -  *  +  En 
I  r>  ria 

where  £n  is  the  energy  of  the  nC^  state  of  the  atom.  The  summation 
in  equation  (2.2)  includes  both  a  sum  over  bound  spin-space  states 
and  an  integration  over  the  continuum  states  of  the  helium  system. 
The  function  Fn(r^t  «*)»  representing  the  nC^  state  expansion 
coefficient,  is  a  onerelectron  scattering  wave-function.  The 
operator  A  anti-symmetrizes  the  product  functions  as  required 
by  the  Pauli  Principle. 

If  the  wave  number  of  the  incident  electron  is  Jf9  , 
in  atomic  units,  then  the  one-electron  scattering  wave-function 
has  the  asymptotic  forms 


j  (2.3) 


F0 (?3  a3)  -  ExpiiAoy)^  +  gx-p  (  ikor3)  a0  (e3)  c r3 


(2.4) 


(2.5) 


where  the  function  denotes  the  plane  wave  incident  along  the 

^  axis  toward  the  target  atom  in  the  zero  bound  state  configuration 
and  an  elastic  outgoing  spherical  wave  in  the  same  spin  state. 

The  spin  function  CTj  denotes  the  familiar  normalized  OC  or 
spin  state  representing  electron  spin-up  and  spin-down,  respectively. 
The  function  Fn(f$C$)  denotes  an  inelastic  spherical  outgoing  wave 
in  the  allowable  spin  states,  Oj  ,  provided  is  real.  If 
is  imaginary,  the  incident  energy  is  less  than  the  energy  required 
for  an  internal  transition  of  the  atom  to  occur,  and  the  boundary 


conditions,  equation  (2.5),  describe  a  decaying  spherical  wave,  the 
sum  of  which  describes  the  polarisation  of  the  initial  atomic  bound 
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state.  The  constraints  on  the  total  system  energy  are  given 
asymptotically  as 


£  =  Ml  ♦  £o 
-  Ml  *  Gn 


(2.6) 


Once  the  detailed  solution  to  the  Schrodinger  equation 
(2.1)  is  available,  the  scattering  amplitudes  }»<•}  can  be  determined 
from  the  asymptotic  boundary  conditions  defined  by  equations  (2.4) 
and  (2.5).  As  shown  in  Mott  and  Massey  (Ref.  19),  the  total  scattering 
cross-section  at  incident  energy  Ml  is  given  as 

Q.(4t2)  =  air  J  f$n(e) I x  SIN  e  dL6  (2.7) 

« 

Close-Coupling  Approximation 

Although  the  expansion  represented  by  equation  (2.2)  is 
formally  exact,  several  practical  difficulties  limit  all  scattering 
calculations.  First,  for  atoms  other  than  hydrogen  or  one-electron 
systems,  the  target  eigen-states  are  not  known  exactly.  Second, 
the  bound  target  states  form  an  infinitely  denumerable  set  of  functions. 
Finally,  the  complete  set  of  target  states  must  include  the  ionization 
continuum.  In  this  case,  the  state  index  n  is  a  continuous  variable 
and  the  number,  though  infinite,  is  not  denumerable. 

In  practice,  only  a  relatively  small  number  of  states 
can  be  included  in  the  expansion  (2.2)  and  are  generally  carefully 
selected.  However,  a  number  of  methods  have  been  developed  to  account 
for  those  terms  neglected  in  the  expansion.  The  basic  approach  of 
these  methods  is  to  replace  the  unperturbed  target  state  wave-functions 
with  pseudo-state  functions  corresponding  to  the  first-order  perturbation 
of  the  target  state  by  a  polarizing  field.  These  pseudo-states 
can  be  constructed  in  terms  of  the  unperturbed  target  wave -functions 
or  other  convenient  basis  functions,  Danburg  and  Karule  (Ref.  21), 

Danburg  and  Geltmann  (Ref.  22).  The  reader  is  referred  to  Mott  and 
Massey  (Ref.  19)  and  Nesbet  (Ref.  23)  for  a  complete  discussion  of 
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the  pseudo-state  method. 

Since  the  purpose  of  this  study  is  the  application  of  a 
variational  principle  to  the  elastic-scattering  of  electrons  from 
helium  in  the  atomic  ground  state,  it  is  convenient  to  adopt  the 
simplest  case  of  the  close-coupling  approximation.  For  incident 
energies  of  the  continuum  electrons  less  than  the  first  excited 
state  of  helium  (19  e.v.),  the  one-state  approximation  is  assumed. 
Thus 


a  %(?,<rt,?x'rA)F0<?3)a3  (2.8) 

where  Yoft* denotes  the  unperturbed  ground  state  helium 
wave -function,  denotes  the  space  part  of  the  incident  electron 

wave-function,  and  <A j  denotes  the  electron  incident  in  the 
arbitrary  spin-up  state. 

Substituting  the  space-spin  hexium  wave- funct ion  and 
anti -symmetrizing  by  a  cyclic  perturbation  of  the  electrons,  we 
have 


»r3^3  )  =  Vof  ri  i0*!  **3  ATa  ,2 

+  %  <3  r,  )F0  (•£>  -  «,  fis 

where  the  spin  function  is  anti-symmetric  with 

respect  to  the  bound  helium  electrons  1  and  2  and  corresponds  to 
a  doublet  total  spin  state  of  the  helium  and  incident  electron. 

is  the  space  part  of  the  helium  wave-function  and  is 
symmetric  with  respect  to  electrons  1  and  2. 

If  we  substitute  the  total  system  wave-function  (2.9) 
into  the  Schrodinger  equation  (2.1),  pre-multiply  by 

sum  over  spin  and  integrate  over 

and  we  obtain 
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Voa  ( ) 


co  oo 


FC(it3)SS  f  ><?">(*.  4  (2.10) 


o  o 


where 


CO  0* 


'oo 


{rS}  =  “If  (-*  -  J  "A  )|M;Cr(^)jx  ^ofr, 
op  r3  r*a  r*a 


and 


Qoo (*)  ?*  ^  )  *  *3)  vf %<Fi fk>  +  **fe(  *0  ^  > 

+  Vn4)^(«4){A  ♦  £ l 

*-  rt  ra.  *3  */a  *7*  *a3  J 

+  (n 


(2.11) 


(2.12) 


In  order  to  obtain  this  result  we  used  Green's  theorem, 
the  symmetry  properties  of  the  electrons  and  the  equation 


OO  OO  r 

J  /  %  (*•  &  I  ^i*+  ^a+— +— "  —  ■*'£o  1  ^o(^ra)o(^  0/4  =  o 

0  0  L  f,  rz  r\x  J 


C2.13) 


This  condition  is  certainly  satisfied  if  is  an 

exact  solution  to  the  helium  equation  (2.3).  Unfortunately,  no 
exact  helium  solution  exists.  However,  equation  (2.13)  can  also  be 
satisfied  if  a  variational  wave-function  for  the  ground  state 
of  helium  is  used  and  the  ground  state  energy  £0  is  given  by 
the  variational  energy.  In  this  study  we  have  adopted  Hylleraas' 
variational  function  given  as 


*o(r,.ra>*  Z*  £xp(-(2rA^Zri)) 

TC 


(2.14) 
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where  Z~  37/lb  and  EgS-AZ9"  and  which  satisfies  equation  (2.13). 

Equation  (2.10)  represents  a  partial  integro-dif ferential 
equation  for  the  scattering  function  Fe(»^)  .  Its  exact  solution, 

which  satisfies  the  boundary  conditions  of  equation  (2.4),  formally 
determines  the  scattering  cross-sections.  For  low-energy  scattering 
it  is  convenient  to  reduce  the  partial  integro-dif ferential  equation 
to  a  set  of  ordinary  integro-dif ferential  equations.  Performing 
a  partial  wave  or  angular  momentum  expansion  on  the  scattering 
function  f^(r)  we  have 

oo 

F0(r)=J.Z  fL(r)  PL  (  COS  ©)  (2.15) 

r  k*o 

where  TlCcos  e)  are  the  Legendre  polynomials.  Substituting 
into  equation  (2.10),  pre-mult iplying  by  ^(Cjos  B)  ,  and  integrating 
over  all  scattering  angles  and  using  the  Hylleraas  helium  wave-function 
we  obtain 

oo  oo 

£-Voo(r3)-ki!z*l))$L<r3)  =J  (fJr.jfJoCr.r^r s)  dr,  clrx  (2.16) 

<*r3  r o  o 

where 


^©o ( *3 )  =  exp  (-a^r3) 


(2.17) 


and 


poo  ( r.  rz  n )  =  »<>Tra  r,  rl r3  %  ( rz ,  r3 )  %  ( **, ,  rx ) 


(2.18) 


and 


VL(r,r%)z  r.Vr,1 


L*l 


-  raL/r,t+' 
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(n<r*) 

(r»<  'i  ) 


(2.19) 


The  boundary  conditions  of  the  radial  function  /fc(r)  are 


tLir)  =  O 

51 N  (M*r  -  Ltt+Il) 

z. 

where  Hl.  is  the  energy-dependent  phase  shift  of  the  L  -partial  wave 
and  is  determined  by  the  exact  solution  to  the  integro-di f ferent ial 
equation  (2.16).  It  can  be  shown,  Mott  and  Massey  (Ref.  19),  that 
the  total  elastic  scattering  cross-section  at  incident  energy  is 


r  •*  o 

r-»  to 


(2.20) 


q(K>)=  2ZZ  (aL+i)  sin*  il 

42  LZ° 

In  the  case  where  the  energy  of  the  incident  electron  is 
small,  only  the  zero  order  phase  shift  is  important  in  the 

calculation  of  the  total  cross-section.  This  was  done  by  Morse 
and  Allis  (Ref.  3),  by  numerical  integration  of  equation  (2.16). 

Setting  L-  o  and  defining 

oo 

K«o  ('•it'y  =  / 
o 


we  have 


(2.23) 


where 


«oo(n.r,)s  l*Z*r,rs  exp(-Zrt)  t*p(-Zr,)  {  (  k?*MZ-l*) 

1  HZ* 


gxp  (-ZZr,)+  txp(-2zrt) 
2*,Z*  zr$z* 


(2.24) 
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III. 


The  Variational  Method 


General  Formulation 


In  order  to  solve  variationally  the  S-wave  integro-differential 
equation,  we  define  the  functional  £  as 


|  (r%)  {  4  Mo  ~  Voo  (ri )  J  $£  (  r* )  of  **3 

«  ae 

-  Jf  fj(n)  K©eCr,»'V(j(r1  o(r3 


(3.1) 


where  is  an  assumed  trial  function. 

Clearly,  if  fo(r)  is  an  exact  solution  to  equation  (3.1), 
fo  (H  ,  then  £>0  .  If  ^o(r)  is  not  exact,  then  jfc.-  constant.  Now 
consider  the  change  in  £  if  is  defined  as 

toM  =  /«(»•)  4 

where  J/o(r)  is  a  first  order  deviation  from  the  exact  solution. 
Substituting,  the  resulting  change  in  is 

°o  a 

^  Sl*°<r3>  +  v'oofr*>jctfo(,3>  cfrf 

3 

-f**<r3>fJfo<r*>koo(r,r-9)dr,clrs  +  #(  l<ff0(r)\*) 

o  o 

where  the  symbol  denotes  integral  terms  of  second  order 

in  the  deviation  from  the  exact  solution. 

Applying  Green's  theorem  to  the  first  terra  in  equation 
(3.3),  we  have 

OB  OO 

I*6<tydL*Jh<v4r9  r  i6f0(rj  <LX  f-a(rs) dlrM 

°  o  dr* 

+  {f0(rJ)^^orrf)  -  ^(r,)  j^f0Cr*>}| 


(3.2) 


(3.3) 


(3.4) 
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where  the  surface  term  is  evaluated  at  the  boundaries.  The  boundary 
conditions  on  &(•*)  are  given  by  equation  (2.20)  for  Ls O  as 

foir'-O 

&(**)  =  A  sin  (  te©r  •+  *lo) 


r-+  o  13.5) 

r  — ►  co 


or  in  different  form 


&(r)«  4,  S/n(A<,r)4aaC^SfJeo»*)  ■*  00  (3.6) 


wl'.ere 


a,  -  A  cos  >io 
A  sin  fio 

and  A  is  a  normalization  constant. 

The  boundary  conditions  on  the  difference  functions  are 

simply 

<ffo(r>  =  o  r-+° 

4fo M  =  4*,  sin  (kc,!*)  +  cos  (k»c )  r-*  co 


(3.7) 


(3.8') 


Substituting  for  the  boundary  conditions,  the  surface  term  reduces  as 


OO 


($><v4<Sfc<rs)  -  4f0(rs)  d  f0<r3)}\  =  +  ~  a,4Qxko 

4*3  «r3  o 


(3.9) 


If  we  impose  on  the  normalisation  constant  A  ,  the  condition 


<W  =  o 


(3.10) 


then  it  follows  that 

+  Qa<fa,  Mo  *  "* 


(3.11) 
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and  the  resultant  change  in  the  functional  &  is 


ML  =  -  *  fdfoty  l  *o  -  v~  <  V }  *©  <ri>  cir3 


(3.12) 


-  IWi'fdfoWKooiWdr'dr*  -f  &(  l *h<r)lZ  ) 
o  o 

Since  the  kernel  function  *oo  (r,r$)  defined  by  equation  (2.24) 
is  symmetric  in  r,  and  fj  we  may  interchange  the  order  of  integration 
and  interchange  the  indices  in  the  third  term  of  equation  (3.12). 

Thus  we  have 


lfo(r3  r  (3.13) 

CO  ^  * 


Substituting  into  equation  (3.12)  we  have 


<«e  =  -  -► 


(3.14) 


where  the  differential  and  integral  terms  vanish  since  £0(r)  is  an 
exac'  solution  to  the  integro-d i f ferent ial  equation  (3.1). 

Since  the  second  order  term  can  be  forced  vanishing  small 
by  choosing  more  and  more  complex  trial  functions  Wr3>  which  better 
approximates  the  exact  solution  f^Cr^)  ,  we  have 

<f(JL+  -  o  (3.i5) 


Thus  the  functional 


is 


stationary  provided 


<W=o. 

Since  the  condition  dfi- O  represents  a  constraint  on 
the  variational  principle  it  can  conveniently  be  incorporated  into 
the  stationary  equation  (3.15)  by  use  of  a  Lagrange  multiplier. 
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From  equation  (3.7)  we  have 


1o  ■  tanT'ca*/**) 


(3.16) 


If  we  require  the  normalization  of  the  asymptotic  solution, 
equation  (3.5)  as  A- 1  ,  which  satisfies  the  condition  6A- O  ,  then 
the  constraint  equation  becomes 


a?  ♦  -  i  *  o 


(3.17) 


Substituting  equation  (3.16)  into  (3.15),  and  introducing 
a  free  parameter  A/A  ,  a  Lagrange  multiplier,  the  stationary  functional 
can  be  written  as 


x-  £.+  -  £(<*?■*** -a;  (3.18) 

Matrix  Variational  Method 

Let  denote  a  trial  function  defined  as 

i9(r)  -  Z  a«4>n(r)  (3.19) 

where  is  an  assumed  set  of  functions,  and  {an}  a  set  of 

undetermined  coefficients.  The  conditions  on  the  function  set  are 


<P±(r)s  sin  (kor) 

<P*(r)  *  cos(Jfer)(i  -  exp  (-a©*-)) 

and 

{<Pn(r)}sO  r  =  O  (ni3) 

{^n(rjJ=o  r-*  oo  (n  >  3) 


where  «o  is  an  assumed  positive  definite  parameter.  We  assume  that 


15 


the  trial  function,  equation  (3.19)  is  an  approximate  solution  to 
the  integro-dif ferential  equation  (2.23).  Clearly  by  construction 
the  trial  function  satisfies  the  boundary  conditions  of  the  exact 
solution,  equation  (3.6). 

If  we  substitute  the  trial  function  into  the  functional 
equation  (3.1),  and  perform  the  differential  and  integral  operations, 
the  functional  £  becomes 


H  N 

X.  :  £  a*  (  Anm  ~  Sum  )  9m 


(3.22) 


where 


A  »»rri  =  J  4n  ( r3 )  {  4  k?  -  Voc  (ra>  \  4m  (r *3 )  alr3 

“  dr.* 


(3.23) 


and 


CD  00 

8nm  s  J  4>n  <r3)  /  4>m  (r, )  Koo  (r,  r3  >  dir,  air 3 
o  o 


(3.24) 


Since  the  functions  {^(^Jare  known,  the  ^  and  B  elements 
are  numbers  corresponding  to  a  given  value  of  the  electron  energy  . 

The  functional  /  at  a  given  electron  energy  reduces  to  a  numerical 
set  of  terms  in  the  undetermined  coefficients  .  Thus  the  functional 

which  is  stationary  under  variations  is 

N  N 

2-  L  L  an(Anm-Q*m)Q.m  ♦  \  (af ♦  a* )  t«N  '(«*)-£  (a, z+a£-i)  0.25) 

mi  m:i  'a.  "*• 


If  we  vary  the  functional  with  respect  to  the  parameter 
set  fM  we  obtain  a  set  of  fi  equations  defined  by 

5  r  O  (3.26) 

6<ln 

Performing  the  differentiation,  the  resulting  set  of  equations 
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in  matrix  form  is 


(6 -**8)0.  -  (A-  Zkolo)™*-  =  O 


where  O  is  an  NxN  syrcmetric  matrix  whose  elements  are 


(3.27) 


8  (  Anm  -  8nm )  ♦  ( 4  Qmn  ) 


(3.28) 


B  is  anW»N  matrix  with  elements  8,^1  ,  B*j  ■ -1  and  the  set 
{•am*®}  ,  Ci  is  the  column  matrix  whose  element  set  is  {a*} 

and  fi\  is  an  NxN  matrix  with  elements  and  {Mum*®}  . 

Formally,  equation  (3.27)  represents  an  NXN  linear  matrix 
problem  which  can  be  solved  for  the  complete  set  of  coefficients,  l*»l  • 
However,  only  the  coefficients  dl,  and  fl^have  physical  significance 
for  determining  the  scattering  phase  shift,  equation  (3.16),  and  thus 
the  collissional  cross-section.  It  is  convenient  therefore  to  reduce 
the  NxM  linear  matrix  problem  to  a  X*X.  eigen-value  problem. 

Let  X  be  a  sub-set  of  the  column  matrix  CL  composed  of 
two  elements  Q,  and  ,  Similarly,  let  Y  be  a  sub-set  column 
composed  of  N-X  element-  .  Thus 


(3.29) 


Defining 

5 


(3.30) 


equation  (3.27)  becomes 

/Cu  Ct»| (X)=u.lX) 
\Ca*  CaaJ\y/  ~\ol 
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where  C„  is  a  3. *3.  array  consisting  of  the  four  elements:  Dm  , 

,  0Zi  • f  he  »  and  .  The  term  C,a  is  a  two  row  by  H  column 

array  containing  the  element  sets{Om}  and  for  (3in<N) 

Similarly  is  an  N  row  by  two  column  array  containing  the  elements 
{°m}  and  or  (3  - ^  —  n)  •  In  particular,  C&i  Is  the  transpose 

of  C*.  Finally,  Caa  is  an  symmetric  array  composed  of 

the  elements  ntm>3. 

Expanding  the  matrix  equation  (3.30)  we  have 


C,(X+  CIXY  =  /<* 

c2/x  +  Cxzy  =  o 


(3.32) 


Operating  on  the  lower  set  of  matrix  equations  by  the  inverse  matrix 
of  the  array,  we  have 


y «  - 


(3.33) 


Substituting  for  Y  we  obtain  a  %*X  eigen-value  problem 


(Cu  -  C*iz  C-Z2  Ca/ )  X  s  X 

Expanding  out  the  Cu  elements  the  eigen-value  problem 

becomes 


(3.34) 


(0,,-£„) 


(  0,2 

(Da»-Hai)  J 


where  we  have  defined  the  value  of  the  £  elements  as 


n*  i,a 

m*  1,2. 


(3.35) 


(3.36) 
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and  denotes  the  am  elei>.  nts  of  the  Cia  C-.1  C,,  matrix  operation. 

The  characteristic  equation  for  the  eigen-values  fA.  is 

simply 

/*-Z  '  /<■  {  (°u  "  j  }  •+  {(  &n  ~  &n  )(  Pa*.  ~  £  ) 

-  (Ola>-£(JL- J^)(  o2f  ~  Ezt+Mo)}  sO 

which  can  be  solved  by  elementary  means. 

The  solution  to  equation  (3.37)  yields  two  eigen-values, 
/•f#and  Mx  since  the  characteristic  equation  is  quadratic  in  Ja.  . 
Similarly  there  exist  two  eigen-vector  solutions,  X,  and  X%  ,  as 
well  as  two  values  of  the  phase  shifts  and  for  each  value 
of  the  incident  energy. 

In  practice  it  is  not  difficult,  in  general,  to  select  which 
values  of  the  phase  shift  correspond  to  a  physical  solution  to  the 
scattering  problem.  This  is  due  to  the  fact  that  the  phase  shifts 
are  defined  as  the  arc  tangents  of  the  ratios  of  the  components  of  the 
eigen-vectors,  equation  (3.16).  Thus  one  set  of  phase  shifts  will 
exhibit  a  uni  form  and  continuous  behavior  as  a  function  ot  energy, 
while  the  second  set  of  phase  shifts  exhibits  a  strongly  divergent 
behavior  in  energy.  Although  both  sets  of  phase  shifts  represent 
the  formal  solution  to  the  quadratic  equation,  the  diverging  set  of 
phase  shifts  can  be  rejected  as  the  physical  solution  to  the 
scattering  problem. 

However,  equation  (3.7)  presents  a  far  more  difficult 
problem  for  specific  regions  of  the  electron  energy.  Since  Ca  ^  Cii  , 
the  roots  of  the  characteristic  equation  can  be  complex.  Since  the 
scattering  phase  shifts  are  real  for  all  values  of  the  electron  energy 
less  than  the  excitation  energy  of  the  helium  atom,  the  variational 
method  will  fail.  In  this  case  both  solutions  to  the  characteristic 
equation  yield  non-physical  results. 


19 


Another  problem  results  due  to  the  quadratic  nature  of  the 
characteristic  equation  and  the  continuous  dependence  of  the  solutions 
on  the  incident  energy,  j£.  While  over  a  broad  range  of  the  energy 
spectrum  only  one  set  of  phase  shifts  represents  the  physical  solution 
to  the  scattering  problem,  for  different  energy  regions  the  physical 
set  can  change.  Thus  we  can  expect  a  small  energy  region  in  which 
both  sets  of  phase  shifts  are  simultaneously  diverging.  In  these 
regions  the  physical  phase  shift  must  be  obtained  by  averaging  the 
two.  It  must  be  emphasised  that  there  is  no  physical  significance 
to  the  anomalous  behavior  of  the  phase  shifts  as  a  function  of 
incident  electron  energies.  In  fact,  the  anomalous  behavior  is  simply 
an  artifact  of  the  variational  method  of  calculation  applied  to 
the  continuum  scattering  problem. 

In  Chapter  IV  we  will  develop  a  computational  algorithm 
for  the  variational  calculation  using  Slater-type  trial  functions 
for  the  trial  function  set  (*>  In  particular  we  will 

also  discuss  an  interpolation  technique  to  obtain  values  for  the  phase 
shifts  at  the  anomalous  energy  values. 
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IV. 


The  Variational  Calculation  Algorithm 


Slater-Type  Functions 

The  first  step  in  the  variational  calculation  is  the 
selection  of  a  specific  set  of  basis  functions  to  represent  the  trial 
solution,  f£(r)  .  Clearly,  the  basis  functions  4>,(r)  and  4>X<r)  , 

defined  by  the  equations  (3.20),  are  constructed  in  order  to  satisfy 
the  boundary  conditions  of  the  exact  solution  to  the  integro-dif ferential 
equation,  equation  (2.23).  However,  the  conditions  imposed  on  the 
basis  set  of  functions  for  are  simply  that  they  vanish  at 

the  boundaries,  equation  (3.21).  A  particularly  simple  set  of 
functions  which  satisfy  the  conditions  are  the  Slater-type  Functions 
defined  as 


<£„(*■)=  r  £xp(-0*r) 


(4.1) 


where  (3 1%  is  a  positive  number  which  is  different  for  each  functional 
term  in  the  expansion,  equation  (3.19).  Thus  the  assumed  form  of 
the  trial  wave-function  is 

N 

f* (r)  =  ax  sin  ( her)  +  a* cos  (k0r)(l- exp (-<V))  +  L*n  rexpi'M-.  2) 

n=3 

Substituting  equation  (4.2)  into  equations  (3.23)  and  (3.24) 
we  can  obtain  explicitly  the  matrix  elements  Anm  and  . 

Appendix  I  is  a  compilation  of  the  f\nm  matrix  elements.  Appendix  II 
is  a  compilation  of  the  Qnm  matrix  elements.  Appendix  III  is  a 
compilation  and  an  evaluation  of  the  definite  integrals  in  which 
the  elements  and  flrtm  are  expressed.  Since  the  evaluation  of 

the  matrix  elements  is  particularly  tedious,  the  choice  of  the 
Slater-Type  Functions  was  made  to  minimumize  the  chance  for  error 
involved  in  the  evaluation.  In  addition,  the  definite  integrals 
resulting  from  the  differentiation  and  integration  of  the  functions 
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have,  in  general,  simple  closed-form  solutions.  Clearly  there  is 

little  restriction  on  the  form  of  the  basis  set  for 

However,  in  order  to  take  advantage  of  the  Variational  Method  as  a 

labor  saving  technique  for  scattering  calculations,  there  exists, 

in  practice,  a  pragmatic  trade-off  in  the  selection  of  the  basis  set. 

Clearly,  with  the  set  of  matrix  elements  given  in  the  appendices, 
we  are  in  a  position  to  generate  a  computational  algorithm  to  calculate 
the  S-wave  phase  shift  for  the  elastic  scattering  of  electrons  from 
helium  in  the  atomic  ground  state. 

Computational  Algorith 

STEP  I:  Initialization  of  Parameters 

There  are  two  types  of  numerical  parameters  needed  in  the 
calculation.  First,  the  values  for  Jfo  and  2  denoting  the  wave-number 
of  the  incident  electron  and  the  ground  state  energy  of  the  helium 
atom  relative  to  the  first  ionization  continuum,  £0  =  ',Z*’  (atomic  units). 

For  the  Hylleraas  wave -function,  equation  (2.14),  the  first  continuum 
is  £0=-a.9M7fc  ,  which  restricts  the  values  of  Jto  between  O  </lt0<  i.tg  . 

Since  we  have  neglected  the  contributions  of  the  excited  states  of 
helium,  the  values  taken  by  must  be  restricted  to  less  than  the 

energy  of  the  first  excited  state  of  the  target  atom.  Thus  we  have 
2*1.697  and  . 

The  second  set  of  parameters  required  are  those  defined  by 
the  basis  set,  fl^and  the  set  CM  ,  depending  on  the  number  of 
Slater-Type  Functions  assumed.  Sinfailam  and  Nesbet  (Ref.  12) 
have  had  particular  success  in  the  more  complex  anomaly-free  Variational 
Method  with  a  decreasing  geometric  progression  of  exponents 
Although  this  selection  is  difficult  to  justify  on  mathematical  grounds, 
we  will  assume  that  the  decreasing  geometric  progression  represents 
a  starting  point  in  the  calculation.  Clearly,  accurate  cross-sections 
will  result  only  after  several  different  runs  are  made  using  different 
distributions  of  the  parameters  (M  .  To  initiate  the  calculation, 
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the  following  set  of  parameters  are  defined  0^sl  ,  i  , 

0*  s  lfin-i)/x  ,  for  o  >3  . 

STEP  II:  Calculation  of  the  Matrix  Elements 

Once  the  numerical  values  Z  ,  Jig  ,  0^  ,  jgg  ,  f&n  are 
specified,  the  matrix  elements  Arttn  and  S/tn\  can  be  calculated  explicitly 
using  the  Appendices  and  the  formulae  for  the  evaluation  of  the 
definite  integrals.  Appendix  III. 


STEP  III:  Construct  the  D&r ,  Matrix 

Once  the  matrix  elements  and Bnajare  known,  we  can  construct 
an  auxiliary  matrix  Dnm  given  by  equation  (3.28).  If  we  set  0‘tx-  D,*“J^ 
and  Oil  =  Da, +4  ,  the  simultaneous  equations  are  defined,  equation  (3,31), 
For  illustrative  purposes  we  will  assume  a  two-free  parameter  trial 
function  which  produces  a  matrix  set.  Thus  equation  (3,31) 

becomes 


D„ 

©ia 

a, 

a. 

d'm 

Da*  0x3  ©a* 

-  A 

o* 

©33  ®3^ 

*3 

Dm 

o*(a  Ohs  &HH 

« 

Oh 

(4.3) 


where 


C,i  - 

(o„ 

Oia  \ 

(4.4) 

1  Dai 

©aa' 

(  D33 

Ok.) 

caa~  j 

l  Dhs 

DqHt 

(4.5) 

s  /  0,s  D,h  \ 

V  Oaj  Dx  / 


and 


^SU  ~ 


IX 


(4.6) 
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In  fact,  the  Cm  matrix  is  independent  of  the  number  of 
basis  functions  used  in  the  expansion  of  the  trial  function  . 


STEP  IV:  Calculate  the  £  Matrix 

At  this  point  the  3.*Z.  eigen-value  problem  is  set  up, 
equation  (3.35).  The  most  difficult  part  of  the  program,  from  a 
numerical  point,  is  the  operation  C/*  C*l  C^a  >  t0  a 

2*Z  matrix  with  elements 


-I  T  , 

Bm  -  Cu  C»a  C ix  I  n 
£ix  *  Ca  Cxx  Ctx\,x 
n  CixCxz  C.IZ  I  ij 

CtxC  CL  *  at  I  xx 


(4.7) 


Clearly  in  the  case  illustrated,  this  only  requires  obtaining 
the  inverse  of  the  matrix,  and  can  be  accomplished  by  elementary 

means,  provided  Cxito  .  For  the  general  N-trial  function  case, 

we  require  the  inverse  of  a  N*2«N-2  matrix.  In  fact,  the  speed 
and  efficiency  of  performing  the  Cxx  matrix  operation  represents 
a  pragmatic  bound  on  the  number  of  basis  functions  used  or  the  range 
of  incident  electron  energy  spanned  in  a  single  machine  run.  Clearly, 
if  OCl  Caar®»  which  is  a  preliminary  requirement  for  the 
operation,  the  Matrix  Variational  Method  completely  fails  to  yield 
a  phase  shift  at  the  specific  value  of  Jcq  .  In  this  case  the  trial 
function  assumed  does  not  represent  the  exact  solution  to  the 
integro-dif ferential  equation  to  first  order  and  the  stationary 
principle  must  he  suspect.  Clearlv,  additional  terms  must  be  added  to 
and  the  particular  energy  region  around  ft*  studied  in  a  later 
machine  run.  In  practice,  if  DPT  Cu  *  O  ,  a  simple  print  of  the 
singular  region  is  only  required.  A  return  to  STEP  II  with  an 
incremental  increase  in  the  value  of  k!o  removes  the  calculation 
from  the  singularity  and  the  program  proceeds. 
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STEP  V:  Solving  the  Eigen-Value  Problem 

All  the  elements  are  now  available  to  solve  the  eigen-value 
problem,  equation  (3,37),  and  which  is  of  the  form 

/lX4  9rjU.  +  C  =  O  (4.8) 

which  has  the  solution 

/ 1  =  -  ir  t/jt-vc  (4.9) 

2. 

•J 

and  yields  the  real  eigen-values  and  ,  proveded  lr  > 

In  this  case  the  eigen-values  are  substituted  back  into  equation 
(3.35)  and  the  real  components  of  the  eigen-vectors  (£,  and 

are  determined.  From  these,  the  needed  phase  shifts  are 
calculated  as  '(^a/Gjland  t<u\  '(4a^a«^a  • 

Figure  1  is  a  schematic  plot  of  the  real  eigen-values  H  * 
and  *lx  as  a  function  of  Jt0  for  electron  scattering  from  hydrogen, 
Shankland  (Ref.  16).  Although  the  basis  set  used  by  Shankland  were 
Gauss-type  functions,  the  general  characteristic  of  the  solutions, 
using  the  Vatiational  Method,  are  independent  of  the  basis  set. 

Thus  Figure  1  illustrates  the  general  nature  of  the  solution  to 
the  eigen-value  problem. 

From  the  figure,  it  is  clear  that  at  a  given  wave-number  , 

only  one  phase  shift  has  converged  to  the  exact  solution  to  the 
scattering  equation.  The  other  phase  shift  is  divergent.  In  practice, 
a  simple  print  and/or  plot  of  the  values  Ht  and  as  a  function  of 
is  .  :i:  to  determine  accurately  i  K  p  ••••Heal  phase  shift 

over  the  entire  energy  spectrum. 

However,  there  exist  anomalous  regions  near  J^~./  and  •  7 
where  neither  solution  or  converged  to  the  exact  solution. 

’  .  j 

Similarly,  there  exists  another  region  where  tr  < 

and  the  roots  of  the  eigen-value  problem  are  Imaginary.  Clearly,  if 


25 


i’xact  Solution 


igure  1.  The  Phase  Shifts  1/  and  'l a  for  Plast  ic  Scatt*. 
from  Hydrogen,  as  a  function  of  Incident 
Wave -Number  . 
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high  accuracy  of  the  phase  shifts  is  not  required  in  the  anomalous 
regions,  then  simple  graphical  or  polynominal  fits  can  be  used  to 
estimate  the  value  of  the  phase  shifts.  If  high  accuracy  is  required, 
a  more  powerful  interpolation  technique  is  needed. 

Interpolation  Technique 

Since  the  general  location  and  nature  of  the  anomalies  are 
known,  a  second  machine  run  is  required  to  study  the  anomaly  in  detail. 

In  this  case  the  range  on  J^e  is  restricted  to  the  anomalous  region. 

Real  Phase  Shifts 

Shankland  (Ref.  16)  proposed  an  empirical  technique  to 
obtain  accurate  phase  shifts  in  the  anomalous  regions,  where 
and  la  are  real,  based  on  the  Stationary  Principle.  Clearly,  the 
entire  development  of  the  Variational  Method  is  based  on  the  fact 
that  the  functional  ,  equation  (3.15),  is  stationary 

with  respect  to  variations  of  the  functional.  Thus  the  value  of 
the  functional  31  has  two  real  values  defined  as 

JL,  s  ( -  10) 

JLXS  (4-n) 

Defining  the  weighed  average  of  the  functional  as 

^  =  ~  Xt  IT a  (4.i2) 

we  have 

i 5  (4-13) 
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However,  sine--  «C-o  we  l'inaLlv  Have 

\  -  *-zn>  -  <• 

X.%- «Ci 

In  order  i  o  ealeu lace  nemer  i  eal  1  y  the  valuer  X,  and  Xj. 
we  require  the  entita  set  o  1.'  coefficients  CLrft;i"-An  >  equation 
(3,22),  For’  unat  ely ,  much  of  the  i;i!oni«i:  ion  required  i already 
calculated  as  part  of  the  e  i  gen-va  lue  problem.  In  tact,  the  (  QN) 

net  of  coefficients  can  be  obtained  by  equation  0.32).  Thus  the 
coefficient  sets  tor  each  >;i\  eu  eigt  m-vector  set  (<j,  aa/( 

and  are  the  solutions  of 

(  JL  ~  C-%3,  C  (  Cji  — yC< i  3L  }  C  ia.  )  Y*  =  O 


(  &-  c'^cJzi  C"-/tzJL)~XC„)Y2:  o 

where  the  matrix  operation  Caa.C(2  is  calculated  as  sub -sect  ion 
of  the  2*2  eigr.n -value  calculation.  Titus  the  values  of  £  and 
can  be  determined. 


1  may’  inary  Plias.e.  Si)  i  Its 

In  this  particular  case,  and  the  roots  o!  the 

t  i gen- value  problem  are  into  tinary,  equation  {',.9)..  fhus  the 
variational  method  fails  to  predict  a  real  phase  shift  in  this 
anomalous  rep, ion. 


A  i  ~/*-i 

conjugate 

shifts  We 
number . 


Since  (lit1  ( wo  roots  are  complex  conjugates  of  each  other. 
,  the  components  of  the  eigen-vectors  will  also  be  complex 
s;  X i  -  X*  •  Tn  order  to  calculat  e  the  complex  phase 
require  a  forma  la  for  the  arctangent  of  an  imaginary 

and  ci  are  real.  then  Abramow i t  h 


Let  <jsx*c^  where  ^ 


2d 


and  Segun  (Kef.  2r>)  snow 


1  -  rvrr  +  J-  tA*  *'  (  a>* 

m 

thus  wo  find  'Ji  -  -  Since 


)  +  4  M  ) 


subs  t. i  Lilt  in);  into 


( * .17) 


t  i  on 


(A  .4)  ,  wo  have 


n(r£AL)  =  Im  <  £^<1 

- - 

which  represents  an  interpolation  formula  for  obtaining  tbo  real 
phase  s-’hift  in  tbo  anomalous:  region  of  the  romp  lex  e  i  gen- ;al  lies , 

In  Figure  2  we  will  outline  the  basic  steps  required  Lo 
compute  the  phase  shift:?  for  elastic  scattering  in  algorithm  term. 


i 


x: 


v. 


Cone lus ion 


The  object  of  this  thesis  was  to  theoretically  formulate 
the  Shankland  Variational  Method  for  application  to  the  electron- 
helium  elastic  scattering  problem.  It  was  found  that  the  scattering 
problem  reduces,  in  the  limit  of  low  electron  incident  energy,  to 
solving  an  integro-dif ferential  equation. 

Instead  of  formally  solving  the  differential  equation,  a 
variational  approach  was  used  to  approximate  the  exact  solution. 

The  results  obtained  by  applying  the  variational  stationary  principle 
were  formulated  into  a  general  matrix  problem.  It  was  shown  that 
the  matrix  problem  redues  to  a  2*2.  matrix  eigen-value  problem. 

An  algorithm  was  developed  to  calculate  the  phase  shifts  for 
the  scattering  electron  at  incident  energy  ,  by  solving  the 

eigen-value  problem.  An  interpolating  technique  was  discussed  which 
can  extract  the  phase  shifts  when  the  eigen-value  problem  gives 
anomalous  results. 

It  appears  from  the  theoretical  formalism  that  the  Matrix 
Variational  Method  provides  a  relatively  simple  technique  for  solving 
complicated  scattering  problems.  The  author  suggests  that  further 
work  be  directed  toward  implementing  the  algorithm  numerically 
and  comparing  the  results  with  known  numerical  solutions  of  the 
integro-dif ferential  equation.  Although  the  physical  limitations  of 
the  cheory  preclude  comparing  the  results  with  experimental  measurements, 
the  author  suggests  that  further  work  be  performed  to  develop  the 
technique  for  polarized  orbitals  and  inelastic  collisions. 
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VPPRh'DI  X  A  i  lio  A  M a t  r :  x  K  lt-tnenl  s 

In  this  appendix  we  list  t  hi-  ft  matrix  elements  def  irsed 
hv  ffjii.'U  >\m  ami  using  t  lie  trial  wave  -  t  une  l  i  <  n  tie  lined  !>v 

equation  (a. 21. 


ft,,  =  ^ziL^zz/ko)  -f  *TLx<,3.z./k») 

ko 

ft,z=  -[< *tX3(«o'ko)  +  AatclL.Cao/ko)*  ^^HU<*o*-^z)/ka)~lLi(^z/J^)j 

•f  1Z{  J3((^>  +  a2)/ie<>)-X3Caz/ii0)jJ 


=  (Pn  +  k£)l ls(0h/Jb>)  ~  TLtifa/Mo)  +  ((fSn  4  AZ)/^  ) 

V  *0  *o 

+  (( fa +  zz  )/*,,) 

4? 


ft*,--  *lL<f(a.Z/k<))-  HJLH((SLZ+«0y*0)  + 

*o 

-  12  3C3  ((.AZ  +  Xol/fa) 

'ko 


/)a^-  - 


^  C 4o/Ao)-]L7  {*«*/**)]  + 

-teo 

-*{  e^(a^«oj_x7^2/^)+3[7((az  +  «^)/k0)} 

+  ^  [  4n  -  JL7  (( 2Z+  Wfa)  +  27  <( az+a<y/Jfc*0?} 

-  12 1 TL7(2Z/Mo  )  -  an 7 +  2r7  (Caz+ac*,)/^ )J 

•fco* 
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/Ian  : 


(M*Ah  {Kd(fin/*o)-  l8  U«e*<3»)/Mo)} 

ko 

+  3/3/)  |  Jtq  (l  fin*  Go)/-ko)  "  JI9  ^  Pn  ^  } 


•+!{  £9^  P*+ a2j/iJo)  -  (in  +  q^^z)  /*<>;  j 

4*0 

+  *zrj[q((fin*2l)/ko)  ~H 9  (( /Sn4 0^  +  52  j/i^} 
■4*o 


Am  =  tilC6((fin+Z Z)/kc)+  *Z  2L5((  $n  +  z7.)/k0) 

*0  k? 


fno. 


*0  -feo 

*  ft||  a  /Sn+^zy^;  -  n9 « /*„ y*,  i] 

i?o 

-  ±  lEqltpn+XQ/Mo)  -X9C(/3n+«(o^^2;/^)} 

■fo» 


/Inm  = 


4^ //3  m  )  _  gffm 

(fin+l 3m)3  (  f3rt  + @>yy)X 

n 

+  _ 

/£n*  /9m  +2Z}3 


_ _ _ 

</9/i+/5«*a  2;x 


/Irm  - 


(pZ+k?)  _  ±_  +  JL  +  Z 
-9/3n3  4/3„  (fa+Z)*  Z(3n*Z)3 . 
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APPRI'TIIX  The  2  Matrix  i.  1  eiucnt  f. 

la  this  appendix  wo  1  i-.t  the  &  matrix  »-i>-ru-nts  defined 
v  equation  (3.24)  and  an  i  nr,  the  trial  wave  -  fane  t  ion  titled  by 
quat inn  (4.2), 


B„=  823jr5u^0)  (4&+«Z-z2)T£*U/.k0)  +  aZ  2L*(3Z/kJ 

Jk%  L 

+  3.TL<>(3Z/fro)  -*  u~2 Z)X(,(2./Me)  +  K3  (a. Z/M0) 

+  8Z3KQ(Z/k<3)X1(  AZ/i^) 

■ko 


Sf2.=  ^3j(^±^^za)Ks(z/^o){lL8(z/ka)-JLa(<z^e)//ie0)l 

^  L  J 
+  ci-Z)Xs(Z/J*0){X'j(z/^)-lL>i((Z->*<>)/*0)$ 

+  (l-Z)JI6(Z/k0)  {  Ke(Z/ko)-Ee(i  2 +  <¥„)//?<,)  } 

+  (Z/ko)H5(Z/ko)  [XqOZ/J^J-KqCCSZ+^V^)  } 

-4-  (z/ko  )  Xs  (3 Z/ko)  {  Eg  <  ZAo>  -  itg  (CZ+%  )/Jko > } 

+  l5iZ/Ao)[JLl  (32/^0)-i:^(C3Z-fc<o)/i«o)[ 

+  lLt,(3zyk0)[x9(z/k0)-Xs{^c(o)/ko)} 

+  {lLd<Z/ko)TL3  (A z/ko)  -  Xsrtz+q^yjfc)  jEgftaz+^'ty} 

-►  {^((Z^)/^)!,  ((az+ct,)/#  -  Ks(z/Mc)JL,  (AZ/ko)} 
+  {X8  (( z-Ho)/ko)X6(AZ+cto)/*to)  -  xa  (2/4%;}  ] 
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Qin  =  ez 3 [ Hs(Vko)  f  (k°  +  hz-zz)  +  <±-z)  + _ zz _ _ 

Mo  Mo  t  (Z  +  fin)*  4fe(3Z+/9»i)* 


(32  /^c) 


l  1  +  xzKs(^/Mq)  (  aJ U(32/ 
JpoOZ+^n)3-!  ^o(2  +  ^n)3  (Z-*-fSn) 


-  a.z gfc(z/^bo)  4  ff5((32*/6*)/Jgo)  +  axt>((zz4p*)'jk>) 

(Z+f!>n)*  Moil*?*)*  (Z +/3*)3 


021=  6Z 


SZ' 3  '(ko  +  HZ-Z*)  X s  ( Z/*o jf  JL 8  ( Z/*o )  - 3C 8  « z+  «o VJeo )  l 

Mi  zMo  1  ; 

+  a-Z  )25  (2/i?o)  {  Xq(Z/Mo)~E3U  Z-kx  0)/ko>] 

+  U-Z)li;6(z^oJ{;Ii:8^z/^>~:IL»((2+OCo>/'ko)  1 

+  ( 2/Jfc)  XiT  ( Z/Mo)  {  3t 8  (3Z/*o> ~T 8  a  3Z+  } 

+  (Z/Mq)  Xs  (3Z /Mo)  {  JL  8  (Z/Mo)  -  %Q  ((  Z  ■+  *0 )  /^)  ] 

+  JLS(ZSMo){T()  (3Z/Mo) -ILtj  ((3Z-t  of^/Mo)  } 

+  1.6  (3Z/Mo){  I g  (Z/ko )  -  JL  9  (( Z+  <*o ) /M0  )  } 

+  Z8  (Z/J?e){l3(aZ/Jte)-  £3  (( *Z +<*<>) /*0)\ 

+  TL  S  U/ko)  [  JL  7  ((zz/Mo)  ~  JL  7  (( 3Z+  «0 ;/*,)} 

-  K*  (z/V  fJL^Uz/Mo)  -  X<j((zz  +  <Xo)/ko)\ 
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&2Z-  8Z3  Ug  +  HZ-Z*)!  (Ze(Z/*o}  -  Zgtfz+cto)/*,))*  } 

ki  L  ***•  1  J 

+  (i-  ax)  { ( Ks ( z/*o) -TLg  (( ) ( JL?  ( Z/ko)  - X?  (( 24  } 

+  (^V{[2L8(Z/feo)  -  I8  «Z+^)/^))(K8<3^o)-IL8«32^)/^))} 

+  a{(!L8(2yjeo)-X8((z^t^)/i^))(lt<,(3z/^)-ir9((3z+of0)/>?<>))} 

+  JLg(Z/kc){jL7(aZ/ko)-^7((^^)/^} 

-  X8«z+°fo>/^{l7  ((azi+aoQ,)/^)^ 

-  iLs(z/k0)fx3(zz/k0)-u:3a  zz+o(0)/jio)} 

+  IE5  ({2+a0>/ko}{'fc3((zz-+cx<>)/ko~TL3(Uz-i- 


£[( 


TEaUAfeo)-  l»((z+^o)l  f  (*£4^z -zaj 


If  (*£-MJZ 
(Z4/3n 


/Sn)3 


+  Cl-  Zj  +  3.Z  -t-  1  1 

(Z*|3n)a  (3Z-fc^J*  (3Z  +  j3n)a' 

+  az[j8  f3z/*o)  -  JLaQz+ebVko)}  -  a  {g9  (z/ka)-J[q  (U+ctpVk,)  \ 

ko  (Z  +  0n)3  (Z+f3n)3 

4  a.|l9(3z/ig0)-I9(3z-<-q0)/^<>?U  f  ITsftaz+ftnMfcO  -ggCCaz^n^o^l 

(Z+£n)a  i?o(Z  +  /3r»;a 


39 


Qni=  gzfosfZ/^f^o  +  ^Z-Z*)  t  U~Z)  t  AZ  1 

^o[  l  (2+^n)3  ^(Z-*-/3r,)a  ^o(3Z+0n)3  ^  (3Z>£r,)a 

+  aZlg(3Z/.Ao)  +  3lU(3ZAfc,>  +  !E8(Z/ito)g3^aZ4^  j^toj 
AolZ-tfir,)^  (Z  +  /3*»)3  i?^* 

+  ag-zJJl6(2/ite)  -  .ZLy  < z/>eo j  +  TLs(2/.k<i)Ie(Uz*f Sn'Ako)  " 
(z+/9/,)3  J?0(z+p„)x  Ve83 


en2=8Z3  \(  ZeiZ'Xo) -JLB«z+«o)/Jto))l  (ieo-f  ^z-z2)  +  d-Z  ) 
r  *»  1  U+/B*;’  (Z-0n)a 


az 


f3Z+«o)3 


+  1  l  -  He(Z/^c)  4  HgfCZt-obJ/ieo) 

(ZZ-t-0«)x)  AfoU+fln)*  Jf0(Z+3n)x 


a.(i~z)£  JLifz/Jto) -IDj  {(z+c(o)/&>)j  +j]Ls(z/Jd0)lLe(C^.x.+jSn)/Ji0)j 
/  -»  .  5  if q  3 


(Z+(Sn)3 


•+  az  f lg  ( 3Z/Jio)  -  IT 8  ((  3Z  »  dpi/ifo )  1 4.  a|'X9(3z/fe)-]I«}(3Z-K<o)/ito)I 
ito  fZ+|9nJ  ^  (Z+Pn)2' 

~^Ifg  ((Z+c<o)/lto)'X8((  3Z+|Sr\+«o?/itojj.  “  £l£>(z/i?o)  *3C5  <UZ+/3„  Vi?o)  j, 

A3 


. 

*f-j  Hs(l  z  +  c(o)/M.o)  •%$  ((az  +  (3n+oro;/ie<>)j 

'Ao  * 
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Bom 


(Z+  /S n)3  { 


a.fi-2) 


(ZL+/3n)3(2  +  /3m)' 


-f  3.^-2)  + _ *2 _ 

(Z+f*n12(Z+(3n)3  (Z*f3n)3(3Z  +  /3r*)3 


-f  ^2.  + _ 5- _ 
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